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ABSTRACT 


Heat  units  are  used  to  characterize  the  interval  between  plant 
and  insect  events  in  terms  of  temperature  unit  accumulations.  Their 
usefulness  is  limited  to  the  temperature  range  in  which  growth  and 
developmental  rates  increase  linearly  with  temperatures  and  do  not 
change  significantly  with  the  age  of  the  organism.  This  publication 
compares  seven  different  methods  of  estimating  daily  heat  units  from 
minimum  and  maximum  temperatures  with  those  calculated  as  the  average 
hourly  heat  units  during  the  day.  The  methods  range  in  complexity 
from  the  simple  mean  temperature  minus  a lower  threshold  to  determining 
the  half-day  areas  that  lie  below  the  minimum-maximum  temperature  line 
by  triangulation  or  sine-curve  by  trigonometry.  At  midtemperature 
ranges  between  the  lower  and  upper  thresholds,  all  the  methods 
obtained  nearly  the  same  number  of  heat  units.  At  lower  temperatures, 
the  simple  mean  method  underestimated  heat  units  and  overestimated 
them  at  higher  temperatures.  The  triangulation  and  sine-curve  methods 
produced  S-shaped  heat-unit  curves  that  tend  to  follow  growth  tempera- 
ture curves.  Microcomputer  algorithms  for  the  seven  methods  were 
listed.  The  triangulation  and  sine-curve  algorithms  were  simplified  to 
provide  maximum  utility  for  predicting  crop  and  insect  developmental 
stages . 


KEYWORDS:  Temperature,  heat  unit,  degree  day,  day  degree, 

triangulation,  sine-curve,  trigonometric  algorithms. 
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INTRODUCTION 


The  growth  and  development  of  plants  and  insects  can  be  characterized  by 
the  number  of  days  between  observable  events,  such  as  cotton  seedling  emergence 
and  first  squares  of  the  duration  of  insect  generations.  The  number  of  days 
between  events,  however,  may  be  misleading  because  growth  rates  vary  with  tem- 
peratures. The  measurement  of  events  can  be  improved  by  expressing  development 
units  in  terms  of  the  temperature  and  time.  The  deviation  between  events  is 
then  based  on  accumulated  degrees  per  unit  time  above  a lower  temperature  re- 
presenting a threshold  of  growth. 

Degree  days  or  day  degrees  can  be  defined  as  days  in  terms  of  degrees  above 
a threshold.  Several  methods  may  be  used  to  calculate  degree  days,  however, 
each  method  obtains  different  values.  Because  of  these  differences,  I will  use 
the  term  heat  unit  (HU)  in  lieu  of  degree  days  for  each  method.  When  using  HU's 
to  time  biological  events,  nonrandom  errors  may  accumulate  because:  (1)  growth 
and  development  rates  of  the  organism  do  not  change  linearly  with  temperature, 
(2)  the  rates  may  vary  with  the  age  of  the  organism  (Wang  I960),* 2  and  (3) 
variables  other  than  temperature  (for  example,  nutrient  status)  may  influence 
the  rates.  In  addition,  (4)  substantial  errors  can  occur  when  lower  and  upper 
threshold  temperatures  are  not  determined  correctly  or  shift  with  the  acclima- 
tization of  the  organism,  (5)  small  inaccuracies  in  temperature  measurements  may 
result  in  substantial  error  when  accumulated  over  an  extended  period,  (6)  the 
calibration  of  growth-stage  HU's  for  actively  transpiring  field  crops  by  using 
temperatures  measured  at  a distance  or  even  in  a standard  instrument  shelter 
near  the  field  may  introduce  large  HU  errors  over  time,  and  (7)  the  experimental 
timing  intervals  that  are  obtained  from  one  method  of  HU  calculations  often 
differ  from  those  obtained  with  other  methods  of  HU  calculations. 


^Plant  physiologist,  Agricultural  Research  Service,  Cotton  Research  Center, 
4207  East  Broadway  Road,  Phoenix,  Ariz.  85040. 

2The  year  in  italic,  when  it  follows  the  author's  name,  refers  to  Litera- 
ture Cited,  p.  14. 
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In  this  bulletin,  I address  the  last  source  of  error  (item  7)  by  comparing 
several  methods  of  estimating  daily  HU's  from  minimum-maximum  temperatures  with 
"observed"  HU’s  in  which  the  hourly  degrees  above  a threshold  temperature  of 
55°F  (12.8°C)  are  summed  from  the  field  observations  and  divided  by  24  hr. 


HISTORY  OF  METHODS 


A critical  review  of  HU's  and  their  use  in  plant  response  studies  was  pre- 
sented by  Wang  (i960).  The  most  commonly  reported  HU  was  the  number  of  degrees 
by  which  the  daily  mean  temperature  exceeds  a lower  threshold  or  base  tempera- 
ture. In  this  method,  the  mean  temperature  was  estimated  as  the  simple  average 
of  the  minimum  and  maximum.  Thus, 


HU  = ((  minimum  + maximum  ) / 2)  - 55 °F 


(1) 


This  is  termed  the  means  method  of  calculating  HU's  (Arnold  1960;  Wang  1960; 
Cross  and  Zuber  1972). 

Lindsey  and  Newman  (1956)  introduced  a method  whereby  HU's  are  calculated 
as  the  triangular  or  trapezoidal  areas  under  the  straight  lines  connecting  daily 
minimum  and  maximum  on  a tempera tur e-time  plot.  Later,  Arnold  (1960)  compared 
the  accuracies  of  calculating  HU's  using  that  triangulation  method  with  the 
means  method,  and  also  a trigonometric  method  whereby  a sine-curve  is  fitted  to 
the  daily  minimum  and  maximum  temperatures.  He  reported  large  errors  in  the 
means  method  and  smaller  errors  in  both  the  triangulation  and  sine-curve 
methods.  Arnold's  comparisons,  however,  involved  only  the  lower  threshold  tem- 
perature for  development,  and  because  only  moderate  temperatures  were  used  he 
did  not  consider  situations  where  high  temperatures  could  retard  development. 


In  cotton  growing  regions,  much  higher  temperatures  are  commonly  recorded 
and  should  be  considered  when  calculating  HU's  (Wang  1960).  Above  an  upper 
threshold  or  optimum  temperature,  the  rate  of  plant  or  insect  growth  may  be  con- 
stant or  may  even  decrease  (Gilmore  and  Rodgers  1958;  Stinner  et  al.  1974; 

Butler  et  al.  1978).  In  HU  algorithms  that  use  temperature  means,  the  maximum 
temperature  is  commonly  limited  to  the  upper  threshold  temperature  (McKinnion  et 
al.  1975;  Gutierrez  et  al.  1975).  In  the  triangulation  (Sevecharian  et  al. 

19-76)  and  sine-curve  methods  (Allen  1976),  however,  the  upper  threshold  has  been 
used  as  an  integral  part  of  the  algorithms. 


REVIEW  OF  METHODS 

Heat  Units  From  Hourly  Temperatures 


Air  temperatures  were  sampled  every  20  min  at  1 m above  the  cotton  crop 
canopy  from  March  to  August  30,  1978,  using  the  data  acquisition  system  de- 
scribed by  Fry  (1978).  Three  readings  were  averaged  hourly  from  which  the  lower 
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threshold  was  substracted.  The  HU's  were  summed  and  divided  by  24  hr.3  The 
HU's  estimated  by  each  of  the  other  methods  are  compared  with  these  daily  ob- 
served HU's. 


Means  Method 


The  means  method  of  calculation  involves  taking  the  mean  of  the  daily  mini- 
mum and  maximum  temperature  and  then  subtracting  a lower  threshold  temperature 
(see  equation  1).  Algorithm  1 in  appendix  A and  figure  1 shows  this  method 
in  which  the  HU's  increase  linearly  with  an  increase  in  the  mean  temperature. 


3Solid-state  devices  that  automatically  sum  HU's  at  10-min  intervals  for 
long  periods  of  time  are  commercially  available.  These  instruments  utilize 
lower  and  upper  threshold  limits. 


MEAN  DAILY  TEMPERATURE  (F) 

40  50  60  70  80  90  100 


MEAN  DAILY  TEMPERATURE  (C) 

Figure  1. — Heat  units  (HU)  calculated  by  the  means,  triangulation,  and  sine- 
curve  methods.  The  continuous  line  represents  the  HU's  estimated  by  the  means 
methods,  algorithm  1 and  2.  The  dashed  line  obtains  HU's  from  the  triangula- 
tion method,  algorithm  4;  and  the  dotted  line,  from  the  sine-curve  method, 
algorithm  5.  The  lower  and  upper  thresholds  were  55°F  (12.8°C)  and  86°F 
(30°C),  respectively,  and  the  swing  was  27°F  (15°C). 


3 


ESTIMATED  HEAT  UNITS  (F) 


Means  Method  With  Upper  Threshold 


A second  means  method  is  listed  in  algorithm  2 in  which  the  maximum  tem- 
peratures are  restricted  to  the  upper  threshold  temperature  (line  220  in  pro- 
gram) before  calculating  the  HU's  as  in  algorithm  1. 


Method  From  Model  SIMC0T  II 


In  the  cotton-plant  model  SIMCOT  II,  McKinnion  et  al.  (1974)  utilized  sepa- 
rate day  (13.2  hr)  and  night  (10.8  hr)  temperature  means  to  calculate  HU’s  as  in 
algorithm  3.  The  combined  HU's  were  similar  to  those  calculated  by  the  means 
method  of  algorithm  2,  but  slightly  shifted  upward  because  of  a greater  number 
of  daylight  hours.  In  addition,  they  used  the  concept  of  the  "physiological" 
day  in  which  a daily  mean  of  79°F  (26°C)  equaled  one  physiological  day. 

Bernhardt  and  Sheppard  (1978)  derived  a similar  physiological  day  from  a mean 
temperature  of  81°F  (27°C)  for  Mexican  beetles. 


Tri angulation  Method 


Sevecharian  et  al.  (1977)  expanded  the  triangulation  method  of  Lindsey  and 
Newman  (1956)  to  involve  an  upper  threshold.  Figure  2 shows  six  possible 
minimum-maximum  temperature  ranges  where  both  thresholds  are  involved  (after 
Allen  1976).  The  shaded  areas  represent  daily  or  half-day  HU's. 

Algorithm  4 calculates  the  areas  according  to  the  flow  chart  in  figure  3. 
In  temperature  ranges  1 and  5 of  figure  2,  the  accumulated  HU's  are  zero  and 
maximum  respectively  (line  530).  When  the  maximum  temperature  is  above  the 
upper  threshold,  the  area  of  the  upper  triangle  D1  (ranges  4 and  6)  is  calcu- 
lated using  the  upper  threshold  H and  the  triangulation  equation  in  line  650. 
Then  the  ranges  2 and  6 that  cross  the  lower  threshold  L are  selected,  and  the 
total  area  is  calculated.  The  area  of  the  upper  triangle  is  then  subtracted 
from  the  total  area  leaving  the  net  HU's  (U).  Range  3 is  calculated  as  in  the 
means  method. 


Sine-Curve  Method  Simplified 


A simplification  of  the  sine-curve  algorithm  is  shown  in  algorithm  5.  The 
range  selection  routine  is  identical  to  that  of  the  triangulation  method 
(algorithm  4).  The  only  difference  between  the  two  algorithms  is  where  the 
trigonometric  equations  replace  the  triangulation  equation  in  the  mathematical 
subroutinue.  These  computer  subprograms  permit  the  choice  of  triangulation  or 
the  sine-curve  algorithms  with  minimal  programming. 

Programmable  hand  calculators  may  be  utilized  to  calculate  heat  units  by 
the  tr iangula tion  and  sine-curve  methods.  Program  listings  for  the  Texas  In- 
strument Company,  Model  TI-59,  are  presented  in  appendix  B. 
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Figure  2. — Temperature  ranges  related  to  a lower  and  upper  threshold. 


TMIN  = DAILY  MINIMUM  TEMPERATURE 
TMAX  = DAILY  MAXIMUM  TEMPERATURE 
L = LOWER  THRESHOLD  TEMPERATURE 

H = UPPER  THRESHOLD  TEMPERATURE 

U = HEAT  UNITS 

M = DAILY  MEAN  TEMPERATURE 

M = ( TMAX+TM I N ) /2 
A = SINE-CURVE  AMPLITUDE 
B,D,  Dl,  T,  Z = TEMPORARY  VARIABLES 

SUBROUTINE  FOR  TR I ANGUL A T I ON 
D = ( TMAX-T ) * *2/ ( 2* ( TMAX-TM IN ) ) 

SUBROUTINE  FOR  SINE-CURVE 
A = TMAX-M  : B = T-M  : Z = ARCS I N ( 8/A ) 
D = ( A*COS(  Z )-B*  ( ( tt/2  )-Z  ) )/tt 


Figure  3. — Flow  chart  for  triangulation  and  sine-curve  algorithms. 
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Sine-Curve  Method  of  Allen 


Arnold  {I960)  introduced  the  concept  of  using  trigonometry  to  calculate 
HU's  from  minimum  and  maximum  temperatures.  Baskerville  and  Emin  (1969)  pre- 
sented the  trigonometric  equations  for  the  sine-curve  treatment.  Allen  {1976) 
published  algorithms  in  the  computer  language  FORTRAN  for  the  six  temperature 
ranges  (fig.  2)  with  and  without  geographical  correction  regressions.  Allen's 
equations  are  shown  in  algorithm  6 with  a geographical  correction  for  Arizona. 


Method  From  Model  KOTTON 

In  the  plant  model  KOTTON  (Gutierrez  et  al.  1975),  the  higher  temperature 
ranges  were  treated  similarly  to  the  means  method  with  an  upper  threshold,  while 
the  lower  temperatures  were  treated  similarly  to  the  sine-curve  method 
(algorithm  7).  One  degree  Fahrenheit  (0.55°C)  was  added  to  each  maximum  tem- 
perature before  proceeding  through  the  algorithm. 


RESULTS 

Theoretical  Differences  Between  Methods 


To  facilitate  a visual  comparison  between  the  two  means,  the  triangulation 
and  the  sine-curve  methods,  HU's  were  calculated  for  each  method  using  a 
constant  diurnal  temperature  swing  of  27°F  (15°C)  and  were  plotted  against  the 
daily  mean  ((min  + max  )/2)  temperature.  The  lower  threshold  was  set  at  55°F 
(12.8°C)  and  the  upper  at  86°F  (30°C)  (fig.  1). 

The  first  means  method  (line  1)  calculates  highest  HU's  at  the  higher 
temperatures  in  range  4.  By  including  an  upper  threshold  as  in  the  second  means 
method  (line  2),  the  HU's  in  range  4 are  lowest.  Both  methods  calculate  the 
same  HU's  in  ranges  1,  2,  and  3. 

The  triangulation  and  sine-curve  methods  (lines  3 and  4)  obtain  symmetrical 
S-shaped  HU  curves  because  the  areas  under  the  minimum-maximum  lines,  or  sine- 
curves  (fig.  2),  change  with  decreasing  rates  as  the  ranges  2 and  4 pass  through 
the  lower  or  upper  thresholds,  respectively.  The  HU's  differ  the  most  between 
the  triangulation  and  the  sine-curve  methods  when  the  mean  temperatures  are  near 
either  the  lower  or  upper  threshold  (table  1).  Figure  1 shows  that  the  triangu- 
lation algorithm  (4)  calculates  fewer  HU's  in  range  2,  and  more  units  in  range  4 
than  does  the  sine-curve  algorithm  (6). 

For  the  other  algorithms  (3,  5,  and  7 ) in  appendix  A,  the  curves  (not 
shown)  tend  to  follow  those  of  the  means  or  sine-curve  methods.  In  the  model 
SIMCOT  II,  the  55  percent  weighting  of  the  day  temperature  over  the  night  (45) 
percent  temperature  caused  the  HU's  to  average  0.1  HU  above  those  calculated  by 
the  means  method  (algorithm  2).  Calculations  by  an  early  KOTTON  model 
(Gutierrez  et  al.  1975)  follow  those  of  the  sine-curve  method  (algorithm  4)  up 
to  the  mean  temperature  at  which  it  switches  to  the  second  means  method. 
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Table  1. — Heat-unit  differences  between  triangulation  and  sine-curve  methods 
as  influenced  by  average  temperatures  and  daily  temperature  swings * 


Average 
tempera  ture 

Daily  temperature  swing 

°F  °C 

20°F(11.1°C) 

30°F(16.7°C)  40°F(22.2°C) 

50°F(27.8°C) 

° F 

°F 

°F 

°F 

35 

(1.7) 

0 

0 

0 

-0.43 

45 

(7.2) 

0 

- .46 

- .93 

-1.35 

55 

(12.8) 

.68 

-1.02 

-1.37 

-1.71 

65 

(18.3) 

0 

- .46 

- .93 

-1.03 

75 

(23.9) 

0 

.36 

.84 

- .85 

85 

(29.4) 

.67 

1.02 

1.36 

1.70 

95 

(35.0) 

.07 

.56 

1.01 

1.42 

05 

(40.6) 

0 

0 

.05 

.53 

iLower  and  upper  thresholds  were  55°F  (12.8°C)  and  86°F  (30°C),  respec- 
tively. 


Table  1 shows  the  calculated  extent  the  HU's  differ  between  the  triangula- 
tion and  sine-curve  methods  when  both  the  mean  temperatures  and  swings  are 
known.  In  addition  to  the  greater  differences  of  HU's  at  the  threshold  tempera- 
tures, the  differences  increased  as  the  daily  temperature  swings  between  the 
minimum  and  maximum  increased. 


Heat  Unit  Conversions 

Mathematical  conversions  of  HU's  between  the  triangulation  and  sine-curve 
methods  were  investigated  but  found  difficult;  however,  some  field  observations 
in  figure  4 were  used  to  indicate  that  there  was  a reasonably  stable  conversion 
factor  (0.97)  for  the  months  of  June  through  September.  (Sine-curve  HU  = tri- 
angulation HU  X 0.97.)  During  these  months,  the  factor  stability  is  obtained 
from  the  relatively  uniform  daily  minimum  and  maximum  temperatures.  Conse- 
quently, daily  changes  in  the  HU's  are  small  when  compared  with  the  changes  that 
occur  during  the  spring  months. 


Observed  Heat  Units  Compared  With  Calculated  Heat  Units 

The  observed  HU's  were  determined  from  daily  sets  of  hourly  field  tempera- 
tures. In  19  sets,  the  temperatures  were  in  range  2 between  January  3 and 
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Figure  4. — Ratio  of  heat  units  calculated  by  the  sine-curve  to  the  triangulation 
methods  for  the  months  of  March  to  September. 


April  8,  1979,  while  50  sets  fitted  range  4 between  April  17  and  October  11, 

1979  and  1980.  Of  the  70  randomly  selected  sets,  one  in  range  6 was  discarded. 
Ranges  1,  3,  and  5 were  not  represented  in  the  sets. 

The  HU's  were  calculated  as  described  in  "Review  of  Methods”  using  thresh- 
old temperatures  of  55°F  (12.8°C)  and  86°F  (30°C).  Rates  of  HU  accumulation  per 
day  were  not  allowed  to  exceed  a maximum  of  the  difference  between  the  thresh- 
olds. Temperature  minimums  and  maximums  were  applied  to  the  different 
algorithms,  as  listed  in  appendix  A with  the  resulting  HU’s  shown  in  table  2. 


Without  an  upper  threshold,  the  first  means  method  overestimated  the  HU’s 
at  higher  temperatures  in  range  4,  and  underestimated  them  at  low  temperatures 
in  range  2.  When  the  maximum  temperatures  were  limited  by  an  upper  threshold, 
as  in  the  second  means  method,  the  HU's  in  range  4 were  underestimated. 


The  routine  in  SIMCOT  II  (algorithm  3)  obtained  nearly  the  same  HU's  as  did 
the  second  means  method.  The  small  differences  in  range  4 results  from  the 
bias  in  favor  of  the  average  day  temperature  for  13.2  hr,  which  was  higher  than 
the  night  temperature  for  10.8  hr. 


HU's  calculated  using  the  triangulation  method  were  not  significantly  dif- 
ferent from  the  observed  HU's  in  both  ranges.  Note  that  in  range  2 the  HU's 
were  slightly  below  the  observed  and  in  range  4 they  were  slightly  above. 
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Table  2. — Mean  daily  heat  units  (HU)  and  t-values  from  observed  and  calculated 
methods  using  hourly  observations  and  minimum-maximum  temperatures * 


Mean 

HU’s 

T-va  lues2 

Method 

Algorithm 

Range 

Range 

Range 

Range 

2 

4 

2 

4 

° F 

°F 

Hourly  observations 

- 

3.3 

13.7 

- 

- 

First  means  w/o  upper  threshold 

1 

2.7 

14.7 

3.81 

6.85 

Second  means  w/upper  threshold 

2 

2.7 

12.2 

3.81 

10.52 

Model  S1MC0T  11 

3 

2.7 

12.3 

3.61 

10.21 

Triangulation 

4 

3.2 

13.8 

1.18* 

1.26* 

Sine-curve  (simplified) 

5 

3.5 

13.4 

2.93 

2.85 

Since-curve  (Allen's  corrections) 

6 

3.8 

13.7 

6.00 

.14* 

Model  KOTTON 

7 

3.5 

12.2 

2.93 

10.52 

Number  of  samples 

19 

50 

iTempera ture  data  were  from  air  sampled  1 m above  the  cotton  crop  at 
Phoenix,  in  1978,  by  the  data  acquistion  system  described  by  Fry  (1978). 
Algorithms  for  calculating  the  HU's  are  described  in  the  text  and  listed  in 
appendix  A. 

2Paired  t-test. 

*Not  significantly  different  from  the  observed  at  P = 0.05. 


The  simplified  sine-curve  method  obtained  HU's  that  were  significantly 
above  the  observed  in  range  2 and  below  in  range  4.  The  above  relationships 
show  that  the  observed  HU's  lie  between  the  HU's  calculated  by  the  triangula- 
tion and  the  sine-curve  methods  in  both  ranges  2 and  4.  The  HU's  from  the  tri- 
angulation method,  however,  are  closer  to  the  observed  as  shown  by  the  lower 
t-va lues . 

Allen's  geographical  corrections  to  the  sine-curve  method  for  Arizona  pro- 
duced an  excellent  fit  of  the  HU's  to  the  observed  in  range  4.  The  same  correc- 
tions, however,  decrease  the  goodness  of  fit  when  applied  to  range  2 as  compared 
with  the  simplified  and  uncorrected  sine-curve  method. 


The  HU's  calculated  by  the  KOTTON  algorithm  were  identical  to  those  from 
the  second  means  method  in  range  4 but  followed  those  from  a sine-curve  method 
in  range  2.  Note  the  large  t-values  from  the  paired  t-test  in  any  ranges  when- 
ever the  means  methods  were  used. 

Of  the  70  original  sets  of  daily  temperatures,  25  sets  were  selected  for 
days  with  100  percent  possible  sunshine.  Plotted  in  figure  5 are  the  observed 
HU's  (continuous  line)  and  the  HU's  for  the  two  means  methods,  the  triangula- 
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MEAN  DAILY  TEMPERATURE  (C) 

Figure  5. — Heat  units  (HU)  observed  and  calculated  from  25  daily  temperature  re- 
cords for  days  with  100  percent  possible  sunshine.  The  continuous  line  repre- 
sents the  observed  HU's  as  summed  from  hourly  field  temperatures.  HU's  esti- 
mated from  minimum-maximum  temperature  are  shown  for  the  means  method  with- 
out an  upper  temperature  limit  for  the  maximum  (x)  and  with  the  upper  limit 
(+),  the  triangulation  method  (A),  and  the  sine-curve  method  (□).  The  average 
temperature  swing  was  28.4°F  (15.8°C)  (SD  = 1.98). 


tion,  and  sine-curve  methods.  HU's  calculated  by  the  triangulation  and  sine- 
curve  methods  approximate  the  observed  HU's  while  the  HU's  calculated  by  the 
means  methods  depart  similarly  to  the  theoretical  relationships  shown  in  figure 

1. 


Data  for  days  with  less  than  100  percent  sunshine  were  excluded  from  the 
plot  because  the  calculated  HU's  departed  from  those  for  sunny  days.  In  range 
2,  clouds  tended  to  interrupt  a smooth  rise  to  and  decline  from  the  maximum  tem- 
perature resulting  in  fewer  HU's  calculated.  In  range  4,  cloudy  and  humid 
nights,  associated  with  cloudy  days,  produced  warmer  nights  and  higher  minimum 
temperatures,  while  the  adjusted  maximum  temperatures  were  constant  at  the  upper 
threshold  with  or  without  the  presence  of  clouds. 
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Comparison  of  Daily  Heat  Units  With  Half-Day  Heat  Units 


To  reduce  errors,  HU's  are  often  calculated  for  half-day  swings  from  mini- 
mum to  maximum  temperatures,  and  then  from  maximum  to  minimum  temperatures 
(Gutierrez  et  al.  1975 ; Allen  1976;  and  Sevecharian  et  al.  1976),  This  treat- 
ment is  most  useful  when  the  daily  mean  temperatures  are  high  and  are  changing 
rapidly;  however,  the  required  number  of  calculations  is  doubled.  Table  3 shows 
the  total  HU's  accumulated  for  the  five  summer  months  using  the  triangulation 
method  in  daily  and  half-day  mode.  The  maximum  accumulated  differences  between 
the  methods  at  any  day  were  not  more  than  5 HU,  and  always  occurred  during  July. 
Daily  positive  errors  tended  to  cancel  the  negative  errors  so  that  near  the  end 
of  the  growing  season  the  accumulated  error  was  small.  For  these  data,  the 
error  for  any  day  was  less  than  1 HU  95  percent  of  the  time. 


Comparison  of  Heat-Unit  Accumulation  Rates 


The  rates  of  HU  accumulations  for  the  months  of  the  year  are  shown  in  table 
4.  The  usefulness  of  the  HU's  for  timing  plant  and  insect  events  during  the 
summer  months  (June  to  September)  may  appear  to  be  no  greater  than  the  use  of 
the  number  of  days.  During  these  months,  the  HU's  are  accumulated  at  a rapid 
rate  with  only  a small  daily  change  of  rate.  If  the  timing  of  events  is  to  be 
compared  between  years,  however,  then  the  use  of  days  alone  may  lead  to  mis- 
interpretation of  the  data  when  some  seasons  are  cooler  and  others  are  warmer 
than  normal.  In  contrast,  during  the  spring  (March-May)  and  autumn  (October- 
November)  months,  the  daily  rates  of  HU  accumulation  are  changing  rapidly.  These 
daily  HU's  are  lower  than  those  accumulated  over  the  summer  months.  Conse- 
quently, the  rapid  daily  changes  show  a high  percentage  of  change  for  these 
periods.  Thus,  the  HU's  furnish  a convenient  timing  aid  of  events  that  are  tem- 
perature-sensitive during  the  spring  and  autumn  months,  and  for  comparisons  of 
results  within  and  between  the  years. 

CONCLUSIONS 

Problems  of  incorrectly  estimating  HU's  can  be  minimized  by  using  the  tri- 
angulation or  the  sine-curve  method  to  calculate  HU's  from  daily  minimum  and 
maximum  temperatures.  These  improved  methods  involve  the  calculation  of  areas 
that  lie  below  either  (1)  a straight  line  or  (2)  a sine  curve  that  extends 
between  the  minimum  and  maximum  temperatures  and  within  the  boundaries  specified 
by  a lower  and  upper  threshold  temperature.  These  areas  represent  HU's  that  tend 
to  follow  the  S-shaped  growth  rate  for  plants  and  insects  more  closely  than  do 
the  linear  HU's  calculated  by  the  means  methods.  Although  the  triangulation  and 
sine-curve  methods  require  more  complex  mathematics  1 calculations  than  the  means 
methods,  look-up  charts,  programs  for  hand  calculators,  and  simple  computer 
programs  are  available  for  both  methods. 
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Table  3.  — Comparison  of  daily  to  half -daily  mode  of  accumulating  heat  units  (HU) 
from  field  temperatures  using  the  triangulation  method"^ 


Total  HU's 

Maximum  difference 

Source 

Period 

Half-daily  Daily 

HU's  Date 

°F 

° F 

°F 

Phoenix  (calc.) 

3/1-8/30 

2114 

2117 

3.8 

7/31/  — 

Tempe  Farm 

4/10-8/30 

2104 

2101 

4.7 

7/12/75 

Phoenix  Airport 

3/2-8/30 

2300 

2302 

2.6 

7/15/78 

Phoenix  Farm 

3/2-8/30 

2141 

2143 

3.3 

7/10/78 

Parker 

3/4-8/30 

1995 

1997 

3.2 

7/29/80 

*Air  temperature  data  were  from  (1)  cubic  polynomial  regression  equation  for 
the  mean  temperature  as  determined  by  the  National  Weather  Service,  Sky  Harbor 
International  Airport,  Phoenix;  (2)  hygr othermograph  in  weather  shelter  at 
Arizona  State  University  Farm,  Tempe;  (3)  monthly  reports  from  daily  weather 
recorded  at  Sky  Harbor  International  Airport;  (4)  air  sampled  1 m above  the 
cotton  crop  canopy  by  the  data  acquisition  system  described  by  Fry  (1968);  and 
(5)  hygr othermograph  in  weather  shelter  at  Bruce  Church  Farm,  Parker,  Ariz. 
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Thble  4. — Mean  daily  temperatures , rates  of  heat  unit  (HU)  accumulation,  and 
mean  daily  rate  of  change  of  accumulations  for  each  month  of  the  year  for 
Phoenix , Ari z . 


Month 

Mean 

tempera  ture 

Mean  rate  of 

HU  accumulation 

Mean  daily 
of  rate 

change 

°C 

°F 

°C 

°F 

°C 

° F 

Percent 

Ja  nuary 

10.7 

51.2 

1.0 

1.8 

0.014 

0.026 

1.47 

February 

12.8 

55.1 

1.9 

3.5 

.041 

.074 

2.14 

March 

15.4 

59.7 

3.4 

6.2 

.076 

.137 

2.21 

April 

19.8 

67.7 

7.1 

12.8 

.162 

.291 

1.71 

May 

24.6 

76.3 

11.6 

20.9 

.129 

.231 

1.11 

June 

29.2 

84.6 

14.9 

26.8 

.040 

.072 

.27 

July 

32.9 

91.2 

16.5 

29.7 

.013 

.023 

.08 

August 

31.7 

89.1 

16.1 

29.0 

- .029 

- .052 

.18 

September 

28.8 

83.8 

14.6 

26.3 

- .035 

- .062 

.24 

October 

22.3 

72.2 

9.5 

17.1 

- .210 

- .377 

2.21 

November 

15.4 

59.8 

3.5 

6.4 

- .121 

- .218 

3.43 

December 

11.4 

52.5 

1.3 

2.3 

- .031 

- .057 

2.50 

Average  temperatures  were  obtained  from  National  Weather  Service,  Sky 
Harbor  International  Airport,  Phoenix.  HU's  were  calculated  according  to  the 
triangulation  method  described  in  the  text,  using  lower  and  upper  thresholds  of 
55°F  (12.8°C)  and  86°F  (30°C),  respectively,  and  a minimum-maximum  temperature 
swing  of  27°F  (15°C). 
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APPENDIX  A:  ALGORITHMS  FOR  METHODS  IN  BASIC  LANGUAGE 

The  programs  on  the  following  pages  are  written  in  simple  BASIC  (with 
multiple  statements),  an  interpreter  language  that  can  be  used  on  most  com- 
puters. Eight  algorithms  are  combined  into  a single  program  for  the  purpose  of 
comparing  the  methods  of  calculation.  Each  algorithm,  however,  will  run  as  an 
individual  routine  requiring  the  minimum  (Tl)  and  maximum  (T2)  temperatures,  and 
the  lower  (L)  and  upper  (H)  thresholds.  Either  temperature  system  (Fahrenheit 
or  Celsius)  may  be  used  after  setting  the  thresholds  to  the  proper  values.  The 
trigonometric  functions  must  be  run  in  the  radian  mode.  If  the  inverse  sine 
function  (ARCSIN)  is  not  available  in  a computer  language,  an  alternative  in- 
verse tangent  (ARCTAN)  may  be  substituted  (see  line  860).  If  necessary,  the 
multiple  statement  line  may  be  broken  into  single  statement  lines  at  each  colon. 
The  SELECT  PRINT  statements  apply  to  the  BASIC  language  for  the  minicomputer 
from  Wang  Laboratories,  Inc.,  Model  2200S,  to  select  output  peripheral  devices. 
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1 SELECT  PRINT  005  : REM  **CRT  DISPLAY** 

2 PRINT 

3 PRINT 

4 PRINT 

5 PRINT 

6 PRINT 

7 PRINT 

8 PRINT 

9 PRINT 
11  PRINT 
12REM 

13  REM  **********  I N I T I HL  I ZE********** 

14  SELECT  PRINT  215<72>  : REM 
Pl=3.  14159  : SELECT  R : REM 


1 1 * * * :+:  * * * * * * * * * * * * * * * * * * * * * * * * :+: :+: :+: :+:  ^ :+; ;+: ;+;  ^ 
HEAT— UNITS  METHODS  COMPARED..  K.  FRV,  7782,  PROG  'HUCALC#  * 
Ll.  S.  D.  A,  AGRICULTURAL  RESEARCH  SERVICE  •+ 

COTTON  RESEARCH  CENTER,  4207  E.  BROADWAY  RD,  PHOENIX,  R2  * 
SOURCES  ARE  IDENTIFIED  IN  SUBROUTINES 
T 1=M  I N I MUM  < TM I N > , T£=MAX  I MUM  < TM  AX  > !J=HE  AT  UN  I TS  < HU  > 

CHANGE  TEMPERATURE  THRESHOLDS  L $<  H IN  LINES  17,  & 19 
* * * * * * * * * * * * :+: * * * * * :+: * * * * * * * * * * * * * * * * * * * * * + * * * * * * * * * :+: * * * * * * * * -■+ 


15 

16 

17 

18 

19 

20 


PRINT 

PRINT 

PRINT 

PRINT 

PRINT 


: PRINT  : PRINT  " 

: PR I NT  " ^THRESHOLDS — LOWER < L > 

" L IN  USE=";L 

" — UPPER  LH!:-  = 30.  O 

" H IN  USE=";H 


♦LINE  PRINTER  OUTPUT* 
♦SELECT  RADIAN  MODE* 
HEAT-UNIT  METHODS  COMPARED" 


= 12. 


i 


55  F 


86  F".i 


: L=55 


: H=86 


21REM 
22  REM 


cL-> 

26 


29 


* * * ;+:  * * R E P E A T RUNS  ST  ART  HERE****** 

PRINT  : SELECT  PRINT  885  : PRINT 
INPUT  "KEY  TMIN",  T1  : INPUT  "KEY  TMAX",  T2 
U=0  : IF  T2>T1  THEN  26  : PRINT  "INPUT  ERROR, 
M=<T1+T2>X2  : REM 
S=T2-T1  : REM 
SELECT  PRINT  215<72> 

PRINT  "TMIN=";  Tl,  "TMAX=";  T2,  " TME AN=  " ; M,  "SWING";  S 
30REM 

31  REM ALGORITHM  1 < 1ST 

108  REM  **D AI L V MEAN  TEMP  WXO  LOWER 
REM  < ADAPTED  FROM  ARNOLD,  I960, 

U=  <!  T 1+T 2 y X2-L  : REM 
IF  U>0  THEN  140  :IJ=8  : GOTO  160 
IF  LKH-L  THEN  160  : REM 
IJ=H-L  : REM 


TMIN>=TMAX"  : GOTO  24 
♦MEAN  TEMPERATURE 
♦TEMPERATURE  SW I NG 


MEANS  METHOD  > 
THRESHOLD**  <FIG  3, 


LINE  1> 


110 
120 
130 
140 
150 
168 
178REM 
188REM 
288  REM 


HQRT.  SCI. 


REM 


76  ; 682  .J 
*HEAT  UNITS 
♦CHECK  FOR  NEGS 
♦CHECK  FOR  >MAX 
♦MAX  HU 


HU 


PRINT  "1  MEANS  WX  LOWER  THRESHOLD",  I NT <U *100+.  5 >7188 


< 2ND  MEANS  METHOD  > 

W7  LOWER  THRESHOLD**  CFIG 


210 

220 

238 

240 

250 

260 

278 


ALGORITHM  2- 

♦+DAILY  MEAN  TEMP 
REM  < ADAPTED  FROM  ARNOLD,  1968,  HORT. 
T5=T2  : IF  T5<=H  THEN  238  : T5=H  : REM 
U= < T1+T5  > 72— L : REM 

IF  U>0  THEN  258  : U=8  : GOTO  270  : REM 
IF  LKH-L  THEN  278  : REM 
U=H-L  : REM 

PRINT  "2  MEANS  W7  LO  & UP  THRESHOLD", 


LINE 


SCI. 


‘6  ; 682 
♦TMAX  UPPER  LIMIT 
♦HEAT  UNITS 
♦CHECK  FOR  NEGS 
♦CHECK  FOR  >MAX  HU 
♦MAX  HU 

I N T < U * 1 8 0 + . 5 > 7 1 8 0 


280REM 
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290REM  RLGORITHM  3 

300  REN  ♦♦FROM  "SI MOOT  1 1 CROP 
310  REM  < ADAPTED  FROM  MCKINION  ET 
320  L5-13.  2 : REM 
330  T5=T 2 : IF  T5<H  THEN  340 
340  M1-T5-T 1 : REM 
350  Lb- L 5/24  : L7=<24-L5>/24  : REM 
360  D5=M1*L6+T1  : N5=M1*L7+T1  : REM 
370  D6=  < D5-L > +L6  : IF  D5>L  THEN  380 
380  N7=<N5-L>+L7  : IF  N5>L  THEN  380 
390  U—D6+N7  : REM 


5 ..  U b D A .•  A R b - S - 5 2 > 

♦DAVLNGTH-  13.  2 HRS 
+TMAX  LIMIT 
♦TEMP  SWING 
♦DRV  & NITE  TIME  FflC 
♦AVG  DAV  & NITE  TEMPS 
: D6=0  : REM  *DAV  HU 
:N7=0  : REM  *NITE  HU 

♦HEAT  UNITS 


MODEL** 
AL.  19 

: T5=H  : REM 


400  PRINT  "3  - S INGOT  1 1 •- - DAY+NITE",  INT<U*100+.  5>/10@ 


410REM 


420REM  ALGORITHM  4 < RECOMMENDED  TR I ANGULATION  METHOD  > 


500 

510 

520 

530 

540 

550 

560 

570 

580 


REM  **TR I ANGULATION  METHOD  SIMPLIFIED,  K. 
REM  < BASED  ON  SEVACHARI AN  ET  AL.  1977.  J. 
M=  < T1+T2  y /2  : REM 

IF  T2<L  THEN  600  : IF  T1>H  THEN  610  : REM 
D1=0  : IF  T2CH  THEN  560  : REM 
T=H  : GOSUB  650  : D1=D  : REM 

IF  T1<L  THEN  570  : U-M-L-Dl  : GOTO  590  : REM 
T=L  : GOSUB  650  : REM 
U=D— D1  : REM 


FRV**  < F I G 3,  LINE  3> 
ECON.  ENT.  70.;  399  > 
♦TEMP  AVG 

♦TEST  FOR  RANGES  1 •&  5 
♦SEP  2 & 3 FROM  4 & 6 
♦AREA  UPPER  TRIANGLE 
♦RANGES  3 & 4 
♦AREA  TOTAL  TRIANGLE 
♦RANGES  2 & 6 


590  IF  U>0  THEN  620  : REM 

600  U=0  : GOTO  620  : REM 

610  U=H-L  : REM 

620  PRINT  "4  TR I ANGULATION 

630REM 


♦CHECK  FOR  NEGS 
♦RANGE  1 

♦MAX  HU,  RANGE  5 
METHOD",  INTCU+100+.  50/100  : GOTO  720 


640  REM  ♦♦TR  I ANGULATION  EQUATION  < SUBROUTINES 

650  D=  < T2-T  > ^2/  i 2*  <.  T2— T1  > : RETURN 

660REM 


670REM ALGORITHM  5 ■:!  RECOMMENDED  SINE-CURVE  METHOD  > 

70O  REM  /♦♦S I NE-CURVE  METHOD  SIMPLIFIED,  K.  FRV++  CFIG.  3,  LINE  4> 
710  REM  t.  ADAPTED  FROM  ALLEN-  1976,  & SEVECHAR I AN  ET  AL.  1977  > 


720  M=  ■:!  T1+T2  > /2  : REM 


+:TEMP  AVG 


730  IF  T2<=L  THEN  800  : IF  T1>H  THEN 

740  D1=0  : IF  T2<H  THEN  760  : REM 

750  T=H  : GOSUB  850  : D1=D  : REM 

760  IF  TKL  THEN  770  : U=M-L-D1  : GOTO 

770  T=L  : GOSUB  850  : REM 

780  U-D-Dl  : REM 

790  IF  U>0  THEN  820  : REM 

800  U=0  .GOTO  820  : REM 

810  IJ-H-L  : REM 

820  PRINT  "5  SINE-CURVE  '^SIMPLIFIED.:' 
830 REM 


810  : REM  ♦TEST  FOR  RANGES  1 & 5 
♦SEP  2 $,  3 FROM  4 & 6 
♦AREA  UPPER  TRIANGLE 
790  : REM  ♦RANGES  3 & 4 

♦AREA  TOTAL  TRIANGLE 
♦RANGES  2 & 6 
♦CHECK  FOR  NEGS 
♦RANGE  1 

♦MAX  HU,  RANGE  5 
",  INTOJ:+:100+.  50/100  : GOTO  920 


840  REM  ♦♦SINE-CURVE  EQUATIONS  '-SUBROUTINE::' ♦♦ 

858  A—T2-M  : B=T-M  : V=B/A  : Z-ARCS I N «:!  V > 

860  REM  ♦ARCS  INC  V>  EQU I V=  Z - A R C T A N < V / S Q R < 1 . 0 0 0 1 - V 1 2 '•  > 

870  D= C A* COS < Z > -B+  < P1/2-Z > > /PI  : RETURN 

880REM 
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S90REM ALGORITHM  6 

900  REM  K**flLLEN'S  SINE-CURVE  W7  CORRECTIONS** 

910  REM  ([CORRECTIONS  FOR  ARIZONA,  ALLEN,  1976,  ENV.  ENTOMOL.  2;  388 > 
920  M=<T1+T2>72  : A=T2-M  :V1=<L-M>/A  : V2=  C K-M  > A 

930  V=C.  08884+A>-.  28187  : REM  +CORREC  FAC  FOR  ARIZONA 

940  IF  T2>L  THEN  960 

950  R=1  : U=0  : GOTO  1060  : REM  GRANGE  1 

960  IF  TKH  THEN  980 

970  R=5  :U=H-L  : GOTO  1060  : REM  GRANGE  5 

980  IF  T2>H  THEN  1020  : IF  T1>=L  THEN  1010 

990  R=2  : F=  < T2-L  > 7 < T2-T1 > : REM  GRANGE  2 

1000  Z1=ARCS  INC  VI  > :Z2=P172  : GOTO  1050 


1010  R=3  :U=M-L+V  : GOTO  I860 
1028  Z2= ARCS  I N < V2  > : IF  TKL 
1038  R=4  : F=  < H-Tl > 7 < T 2-T 1 > : 
1048  R=6  : F=  < H-L  > 7 < T2-T1  > :Z 
1050  U=  < < < M-L  > •+■  < Z2-Z1  > + C COS  C 
I860  PRINT  "6  SINE-CURVE  <W7 


: REM  GRANGE  3 

THEN  1048 

Zl=-Pl/2  : GOTO  1050  : REM  GRANGE  4 
1= ARCS  I N <:  VI  > : REM  GRANGE  6 

Z1 > -COS C Z2 > > +A+ < H-L > * < P1/2-Z2  > > +17P1 > + < V+F > 
ALLEN - S COREO  ",  INTCU+180+.  5 >7100, 


1870  PRINT  "TEMP  RANGE=".;  R 


1080 REM 

1090REM ALGORITHM  7 


1180  REM  Z++FROM  'KOTTON'  CROP  MODEL** 

1118  REM  < ADAPTED  FROM  GUTIERREZ  ET  AL.  1975.  J.  ENV.  ENT.  4.;  125 > 
1120  F5=l  : REM  ^CONVERT  FACCF5>=1.  8 C,  1.  0 F 

1130  T5=T2+F5  : IF  T5CH  THEN  1148  : T5=H 
1148  U=0  : IF  T5-10L  THEN  1188 

1150  T=THT5-C2*L>  : IF  TKL  THEN  1160  : U=T72  : GOTO  1188 
1168  V=  < T7 < T1-T5 > > : Z=ARCS I N < V > 

1178  U= < 17 C 2*P1 > * < < T5-T1 > *COS C Z > -T+Z > > + < T74 > 

1180  PRINT  "7  ' KOTTGN  ' < MEANS7S I NE  COMBO >",  I NT  ■:  0+188+.  5 >7188 
1190REM 

1208  GOTO  22  : REM  ^REPEAT  RUN 
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HEAT-UNIT  METHODS  COMPARED 


♦THRESHOLDS — LOWER (L*  = 12.  8 C 

— UPPER ( H * — SB.  U U 


55  F L IN  IJSE=  55 
86  F H IN  USE=  86 


TM I N=  40  TMAX*  70 

1 MEANS  14/  LOWER  THRESHOLD 

2 MEANS  W / LO  & UP  THRESHOLD 

3 "SI MOOT  II'-  DAV+NITE 

4 TR I ANGULATION  METHOD 

5 SINE-CURVE  (SIMPLIFIED* 

6 SINE-CURVE  (W/  ALLEN'S  COREO 

7 'KOTTON'< MEANS/SINE  COMBO* 


TMEAN=  55 
0 
0 


4.  77 

5.  3 
5.  19 


SWING  30 


TEMP  RANGE=  2 


TM  I N=  55.  5 TMAX*  85.  5 

1 MEANS  W/  LOWER  THRESHOLD 

2 MEANS  14/  LO  & UP  THRESHOLD 

3 'SI MOOT  II'-  DAV+NITE 

4 TR I ANGULATION  METHOD 

5 SINE-CURVE  (SIMPLIFIED* 

6 SINE-CURVE  <14/  ALLEN'S  COREO 

7 ' KOTTON  " < MEANS/S I NE  COMBO * 


TMEAN*  70.  5 
15.  5 
15.  5 
15  65 
15.  5 

15.  5 

16.  55 
15.  75 


SWING  38 


TEMP  RANGE = 3 


TMIN=  71  TMAX*  101 

1 MEANS  14/  LOWER  THRESHOLD 

2 MEANS  14/  LO  & UP  THRESHOLD 

3 ''SI MOOT  II'-  DAV+NITE 

4 TR I ANGULATION  METHOD 

5 SINE-CURVE  (SIMPLIFIED* 

6 SINE-CURVE  (W/  ALLEN'S  COREO 

7 ' KOTTON ' ( MEANS/S I NE  COMBO  * 


TMEAN*  86 
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SWING  30 


TEMP  RANGE*  4 


TM I N=  87  TMAX=  112 

1 MEANS  14/  LOWER  THRESHOLD 

2 MEANS  14/  LO  & UP  THRESHOLD 

3 'SI MOOT  II'-  DAV+NITE 

4 TR I ANGULATION  METHOD 

5 SINE-CURVE  (SIMPLIFIED* 

6 SINE-CURVE  <14/  ALLEN'S  COREO 

7 ' KOTTON  '•  ( MEANS/S  I NE  COMBO  * 


TMEAN*  99.  5 
31 
31 

31.  5 
31 
31 
31 


SWING  25 


TEMP  RANGE*  5 


TMIN=  50  TMAX*  90 

1 MEANS  14/  LOWER  THRESHOLD 

2 MEANS  14/  LO  A UP  THRESHOLD 

3 'SI MOOT  II'-  DAV+NITE 

4 TR I ANGULATION  METHOD 

5 SINE-CURVE  (SIMPLIFIED* 

6 SINE-CURVE  <14/  ALLEN'S  COREC* 

7 ' KOTTON  ' (MEANS/SINE  COMBO* 


TMEAN*  70 
15 
13 

13.  18 
15.  11 

15.  22 

16.  38 
13.  8 


SWING  48 


TEMP  RANGE*  6 
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APPENDIX  B:  HEAT-UNIT  PROGRAMS  FOR  POCKET  CALCULATOR 


HU’s  may  be  calculated  and  accumulated  on  the  programmable  pocket  calcula- 
tor, Texas  Instrument  TI-59 , with  the  following  program.  Either  the  triangula- 
tion or  the  sine-curve  subroutine  may  be  attached  to  the  main  program,  which 
selects  the  proper  temperature  range  as  shown  in  figure  2.  Print  commands  are 
included  in  the  listing  for  use  with  a printer  but  are  not  necessary  for  opera- 
tion. With  the  aid  of  the  HU  flow  chart  in  figure  3,  this  program  may  be 
adapted  to  other  programmable  calculators. 
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TITLE  HEAT-UNIT  METHODS PAGE  1 OF  3 Tl  PfOgfOmmable 

PROGRAMMER  K.  FRY,  USDA,  ARS  qaTEAUGUST  1982  PfOgfQm  R@COfd 

Partitioning  (Op  17)  1417,9.5,91  Library  Module  NOT  USED Printer  PC-  1 00A  Cards  1 SIDE 


USER  INSTRUCTIONS 


STEP 

PROCEDURE 

ENTER 

PRESS 

DISPLAY 

1 

LOAD  PROGRAM 

2 

INITIALIZE  (HU  SUM  = 0) 

2ND 

E' 

THUP 

3 

KEY  MINIMUM  TEMPERATURE 

TMIN 

A 

TMIN 

4 

KEY  MAXIMUM  TEMPERATURE 

TMAX 

B 

TMAX 

5 

HEAT  UNITS  - MEANS  METHOD 

C 

HU 

6 

PRINT  AVERAGE  TEMPERATURE 

2ND 

A ' 

TAVG 

(STEPS  5 & 6 MAY  BE  SKIPPED) 

7 

HEAT  UNITS  - TR I ANGULA T I ON 

D 

HU 

(SINE-CURVE  METHOD  MAY  BE  SUBSTITUTi 

:d 

BY  CHANGING  SUBROUTINE  ADDRESSES  AT 

LINES  72-4  AND  99-101  TO  SBR  196) 

8 

PRINT  ACCUMULATED  HEAT  UNITS 

E 

HU  SUM 

(THE  PROGRAM  WILL  OPERATE  WITH  OTHEI 

t 

LOWER  AND  UPPER  TEMPERATURE  THRESHOt 

.DS 

BY  CHANGING  THE  VALUES  IN  LINES  146- 

■149 

AND  1 56-  159). 

USER  DEFINED  KEYS 

DATA  REGISTERS  ( f.M)  US  ) 

LABELS  (Op  08) 

a TMIN 

o THUP-THLO 

0 

[INV] Qnx] [CE]  [CLR]  (*Lt]  _ 

b TMAX 

i TMIN 

1 

|7t]  [Vx]  [sfo]  [roI  |sum]  IT1! 

c HU  (MEANS) 

2 TMAX 

12  HU  SUM  (MEANS 

nn  _ m _ m _ m _ umi  _ m _ 

D HU  (TRIANG) 

3 TAVG 

13  HU  SUM  ( TR I AN< 

- pn  rn  irst]  no  __  ists)  _ rn  _ 

E HU  SUM 

4 D 2 

4 

®_c=]  w na  n ei 

■Bl  UTl  IB1  MM  WHS  ITS 

a AVERAGE  TEMP 

5 SINE-CURVE  EQ 

5 

Mffl  iixl  Ihi  ID  IK  H 

B' 

6 

6 

121  ES  K3  8SH  El  111 

C' 

7 " 

7 

arcn  si  gg  ga  ng  mu 

D' 

8 LOWER  THRESH. 

8 

in  m ei  mi  a ii 

e INITIALIZE 

9 UPPER  THRESH . 

9 

in  sa 

FLAGS  0123456789 


_ PROGRAM  DESCRIPTION  _ 

THIS  PROGRAM  WILL  CALCULATE  DAILY  HEAT  UNITS  BASED  ON  MINIMUM  AND  MAXIMUM 
TEMPERATURES  USING  THE  MEANS,  TR I ANGUL A T I ON , OR  THE  SINE-CURVE  METHOD. 

_ I NE  S 0-11  INPUT  THE  TEMPERATURE  DATA,  LINES  12-33  CALCULATE  HEAT  UNITS  BY 
THE  MEANS  METHOD,  LINES  34-120  SELECT  TEMPERATURE  RANGES  AND  CALCULATE  HEA 
JNITS  USING  EITHER  THE  TR I ANGUL A T I ON  (SUBROUTINE  LINES  178-195)  OR  THE 
SINE-CURVE  METHOD  (SUBROUTINE  LINES  196-239).  INITIALIZATION  AND  PRINT-OU 
UTILITIES  ARE  LISTED  FROM  LINE  121  TO  167. 
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title — heat-unit  methods page _2 — of  — 3 — j|  Programmable  firn 

PROHRAMMFR  K.  FRY,  USDA,  ARS DATE Coding  Form  V 


LOC 

CODE 

KEY 

COMMENTS 

LOC 

CODE 

KEY 

COMMENTS 

LOC 

CODE 

KEY 

COMMENTS 

0 

76 

LBL 

KEY 

0 1 

0 1 

TMAX  ? 

1 10 

25 

CLR 

RANGE  1 

1 1 

A 

MINIMUM 

1 0 

1 0 

GOTO  110 

1 1 1 

AA 

SUM 

OUT 

A 2 

STO 

TEMP 

A3 

RCL 

13 

1 3 

0 1 

01 

09 

09 

99 

PRT 

99 

PRT 

( TMIN ) 

32 

X 1 T 

98 

ADV 

9 1 

R/S 

60 

A3 

RCL 

9 1 

R/S 

STOP 

76 

LBL 

KEY 

0 1 

01 

1 1 6 

A3 

RCL 

RANGE  5 

1 2 

B 

MAXIMUM 

77 

GE 

TMIN  >= 

00 

00 

A 2 

STO 

TEMP 

0 1 

0 1 

THUP  ? 

6 1 

GTO 

02 

02 

1 6 

16 

GOTO  116 

01 

0 1 

1 0 

99 

PRT 

( TM AX ) 

A3 

RCL 

120 

1 1 

1 1 

9 1 

R/S 

02 

02 

121 

76 

LBL 

PRINT 

76 

LBL 

MEANS 

22 

INV 

16 

A ' 

TAVG 

1 3 

C 

METHOD 

77 

GE 

TMAX  < 

7 1 

SBR 

A3 

RCL 

W/O  UPPEF 

00 

00 

THUP  ? 

01 

01 

00 

00 

THRESHOLC 

70 

77 

77 

GOTO  77 

5 

68 

68 

32 

X £ T 

32 

X 

AV 

H 

99 

PRT 

7 1 

SBR 

GET 

72 

7 1 

SBR 

GET  AREA 

98 

ADV 

0 1 

01 

( TAVG ) 

0 1 

01 

OF  UPPER 

9 1 

R/S 

68 

68 

7 A 

78 

78 

TRIANGLE 

129 

76 

LBL 

PR  I NT 

20 

75 

- 

A 2 

STO 

( D 1 ) 

130 

1 5 

E 

HU  SUMS 

A3 

RCL 

OA 

OA 

A3 

RCL 

08 

08 

77 

A3 

RCL 

12 

12 

95 

= 

08 

08 

99 

PRT 

22 

INV 

MAX  HU 

32 

X ^ T 

A3 

RCL 

77 

GE 

LIMITED 

80 

A3 

RCL 

13 

1 3 

00 

00 

TO 

0 1 

01 

99 

PRT 

29 

29 

H-L 

22 

INV 

98 

ADV 

32 

X £ T 

77 

GE 

TMIN  < 

Q 1 

» /s 

29 

AA 

SUM 

00 

00 

THLO  ? 

1 39 

76 

LBL 

INITIAL- 

30 

1 2 

12 

97 

97 

GOTO  97 

1 AO 

10 

E' 

I ZE 

99 

PRT 

A3 

RCL 

25 

CLR 

98 

ADV 

03 

03 

A 2 

STO 

9 1 

R/S 

75 

- 

1 2 

12 

76 

LBL 

A & 

32 

X 

AV 

H 

A 2 

STO 

1 A 

D 

AREA 

90 

75 

_ 

1 3 

13 

7 1 

SBR 

METHODS 

A3 

RCL 

1 A6 

05 

5 

LOWER 

0 1 

0 1 

OA 

OA 

05 

5 

THRESH- 

68 

68 

GET  TAVG 

95 

rs 

HU  = 

93 

OLD 

A 2 

STO 

6 1 

GTO 

TAVG  - 

1 A9 

00 

0 

° F 

AO 

03 

03 

0 1 

0 1 

THLO  - 

150 

A 2 

STO 

25 

CLR 

06 

06 

DL 

08 

08 

A 2 

STO 

97 

A3 

RCL 

99 

PRT 

OA 

OA 

08 

08 

9 A 

+ / - 

A3 

RCL 

99 

7 1 

SBR 

GET  AREA 

A 2 

STO 

02 

02 

100 

0 1 

01 

OF  TOTAL 

00 

00 

32 

X S T 

101 

78 

78 

TRIANGLE 

1 56 

08 

8 

UPPER 

A3 

RCL 

75 

_ 

( D ) 

06 

6 

THRESH- 

0 1 

01 

A3 

RCL 

93 

. 

OLD 

50 

77 

0 1 

10 

A3 

08 

77 

GE 

0 1 

10 

RCL 

08 

GE 

TMIN  > = 

TMAX  ? 

GOTO  110 

THLO  > = 

1 06 

OA 

95 

29 

77 

0 1 

1 1 

OA 

HU  - D-D  1 

(CLEAR  T) 

HU  > 0 7 

GOTO  111 

1 59 

00 

0 

°F 

CP 

GE 

0 1 

1 1 

MERGED  CODES 

62  m mg  72  [sToj  m 83|gtoi  mi 

63  ma  mi  73  @ mi  84  0 m 

64  mi  mi  74  no§  mi  92  [By]  ^ 

Texas  Instruments 

1 NCOR  POR  AT  E D 

‘977  Texas  Instruments  Incorporated 
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TITLE  heat-unit  methods 


PAGE  _3 OF 3 

PROGRAMMER  K ■ FRY»  USDA»  ARS DATE 


Tl  Programmable 

Coding  Form 


LOC 

CODE 

KEY 

COMMENTS 

LOC 

CODE 

KEY 

COMMENTS 

LOC 

CODE 

KEY 

COMMENTS 

160 

42 

STO 

42 

STO 

09 

09 

07 

07 

99 

PR  T 

43 

RCL 

44 

SUM 

THUP  - 

06 

06 

00 

00 

THLO 

65 

X 

98 

ADV 

220 

43 

RCL 

70 

RAD 

07 

07 

1 $7 

?l 

R/S 

39 

COS 

168 

43 

RCL 

AVERAGE 

75 

— 

0 1 

0 1 

TEMP 

53 

( 

170 

85 

+ 

SUBR  . 

43 

RCL 

0 

43 

RCL 

05 

05 

02 

02 

65 

X 

2 

95 

= 

53 

( 

55 

89 

TT 

5 

02 

2 

230 

55 

5 

6 

95 

— 

02 

2 

8 

92 

RTN 

75 

- 

7 

178 

94 

+ / - 

TRI ANGU- 

43 

RCL 

9 

85 

+ 

L AT  I ON 

07 

07 

9 

180 

43 

RCL 

SUBR  . 

95 

= 

o 

02 

02 

55 

4 

95 

— 

89 

TT 

33 

X2 

9 5 

_ 

4 

55 

2 39 

92 

RTN 

5 

53 

( 

02 

2 

65 

X 

53 

( 

43 

RCL 

1 90 

02 

02 

75 

_ 

2 

43 

RCL 

0 1 

0 1 

95 

= 

195 

92 

RTN 

196 

75 

- 

SINE- 

43 

RCL 

CURVE 

03 

03 

SUBR  . 

95 

= 

200 

42 

STO 

05 

0.5 

55 

4 

53 

( 

43 

RCL 

02 

02 

75 

- 

43 

RCL 

03 

03 

54 

) 

2 1 0 

42 

STO 

MERGED  CODES 

62  KU  ua 

72  [stoi  mn 

83  [gto]  Hill 

63  m HI 

73  [rclI  US 

84o  m 

9 5 

= f 

64 12|  gjg 

74  [SUM]  Bn 

92[in«|  [sbw! 

2 2 

I N V 

Texas  Instruments 

38 

SIN 

INCORPORATED 
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